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Abstract  Manifestations  of  the  11 -year  solar  cycle  and  longer  time-scale  variability  in  the 
heliosphere  and  cosmic  rays  are  considered.  We  briefly  review  the  cyclic  variability  of  such 
heliospheric  parameters  as  solar  wind  speed  and  density  and  heliospheric  magnetic  field, 
open  magnetic  flux  and  latitude  variations  of  the  heliospheric  current  sheet.  It  is  discussed 
whether  the  local  in-situ  observation  near  Earth  can  represent  the  global  3D  heliospheric 
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pattern.  Variability  of  cosmic  rays  near  Earth  provides  an  indirect  useful  tool  to  study  the 
heliosphere.  We  discuss  details  of  the  heliospheric  modulation  of  galactic  cosmic  rays,  as 
recorded  at  and  near  Earth,  and  their  relation  to  the  heliospheric  conditions  in  the  outer 
heliosphere.  On  the  other  hand,  solar  energetic  particles  can  serve  as  probes  for  explosive 
phenomena  on  the  Sun  and  conditions  in  the  corona  and  inner  heliosphere.  The  occurrence 
of  major  solar  proton  events  depicts  an  overall  tendency  to  follow  the  solar  cycle  but  in¬ 
dividual  events  may  appear  at  different  phases  of  the  solar  cycle,  as  defined  by  various 
factors.  The  solar  cycle  in  the  heliosphere  and  cosmic  rays  depicts  a  complex  pattern  which 
includes  different  processes  and  cannot  be  described  by  a  simple  correlation  with  sunspot 
number. 

Keywords  Heliosphere  •  Cosmic  rays  •  Solar  energetic  particles  •  Solar  activity 

List  of  Abbreviations 


AU 

Astronomical  unit  (the  mean  Sun-Earth  distance,  ^1.5  x 

CME 

Coronal  mass  ejection 

CR 

Cosmic  rays 

GEE 

Ground-level  enhancement  of  cosmic  rays 

GCR 

Galactic  cosmic  rays 

HCS 

Heliospheric  current  sheet 

HME 

Heliospheric  magnetic  field 

IC 

Ionization  chamber 

MHD 

Magnetohydrodynamics 

NM 

Neutron  monitor 

OME 

Open  magnetic  flux 

SEP 

Solar  energetic  particles 

SPE 

Solar  proton  events 

SSN 

Sunspot  number 

1  Introduction 

The  heliosphere  is  a  cavity  in  the  interstellar  plasma  controlled  by  the  solar  wind  contin¬ 
uously  emitted  by  the  Sun  and  by  heliospheric  magnetic  field  (HME)  ultimately  produced 
by  the  solar  dynamo  (see,  e.g,  recent  reviews  by  (Owens  and  Eorsyth  2013;  Balogh  and 
Erdos  2013).  The  exact  size  of  the  heliosphere  may  vary  in  the  course  of  solar  cycle  (e.g., 
Pogorelov  et  al.  2013,  and  references  therein).  It  consists  of  three  main  regions  separated  by 
clear  boundaries. 

-  The  termination  shock  bounds  the  region  of  supersonic  solar  wind  with  frozen-in  HME. 
This  is  an  approximately  spherical  region  with  a  radius  of  about  90  AU.  The  helio¬ 
sphere  inside  the  termination  shock  does  not  “know”  (in  the  MHD  sense)  about  the 
presence  of  anything  outside.  The  termination  shock  was  discovered  by  Voyager- 1  at 
the  distance  of  94  AU  in  2004,  and  by  Voyager-2  at  84  AU  in  2007  (Stone  et  al.  2005, 
2008). 

-  The  region  beyond  the  termination  shock  is  still  filled  by  the  solar  wind  and  HME,  but 
those  are  distorted  by  the  presence  of  the  interstellar  wind.  This  region  is  bounded  by  the 
heliopause  which  is  non-axisymmetric.  Its  size  is  about  122  AU  in  the  nose  direction,  as 
discovered  by  Voyager- 1  in  2012  (Stone  et  al.  2013).  The  size  in  the  tail  direction  is  not 
precisely  known. 
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-  Interstellar  wind  is  somewhat  disturbed  beyond  the  heliopause  forming  an  interface  re¬ 
gion,  whose  exact  size  is  not  well-known.  An  existence  of  a  bow  shock  as  an  inter¬ 
face  between  disturbed  and  undisturbed  interstellar  wind  was  proposed,  but  recent  data 
disputes  that  (McComas  et  al.  2012;  Zieger  et  al.  2013).  Beyond  this  boundary,  inter¬ 
stellar  wind  does  not  “know”  about  the  presence  of  an  obstacle  caused  by  the  solar 
wind. 

The  heliospheric  structure  and  conditions  can  be  probed  by  not  only  dedicated  spacecraft 
but  also  by  cosmic  rays  traversing  the  heliosphere.  Cosmic  rays  near  Earth  are  of  two  main 
types:  energetic  (10^-10^®  eV/nucleon)  galactic  cosmic  ray  (GCR)  particles  existing  all  the 
time  but  modulated  by  the  solar  magnetic  activity  within  the  heliosphere;  and  less  energetic 
(10^-10^®  eV/nucleon)  solar  energetic  particles  (SEP)  which  occur  as  sporadic  solar  proton 
events  (SPEs),  when  the  flux  of  such  particles  can  be  enhanced  by  many  orders  of  magnitude. 

Since  the  solar  magnetic  activity  depicts  a  dominant  11 -year  cycle  driven  by  the  solar 
dynamo  (Hathaway  2010),  the  heliosphere  also  changes  with  the  11-year  solar  cycles.  Ac¬ 
cordingly,  the  flux  of  energetic  particles  near  Earth  also  varies  both  within  the  1 1 -year  solar 
cycle  and  on  the  long-term  scale.  Here  we  present  a  brief  review  of  the  solar  cycle  as  re¬ 
flected  in  different  aspects  of  the  heliospheric  and  cosmic  ray  variability  near  the  Earth. 
Other  relevant  reviews  can  be  found,  e.g.,  in  Balogh  et  al.  (2008). 

Although  this  paper  is  devoted  to  the  near-Earth  observations,  it  is  necessary  to  men¬ 
tion  invaluable  contributions  of  distant  space  missions  to  understanding  of  the  heliosphere 
3D  structure  and  dynamics.  Pioneer  10  (1972-2003),  Pioneer  11  (1973-1995)  reaching 
^80  AU,  and  Voyager  1  and  2  (both  launched  in  1977  and  now  exploring  the  outer  helio¬ 
sphere  and  beyond)  provide  unique  data  exploring  the  extreme  far  regions  of  the  heliosphere. 
The  space  probe  Ulysses  (1990-2009)  explored  space  outside  the  ecliptic  plane,  performing 
three  scans  over  the  solar  polar  regions  in  1994-1995,  2000-2001,  and  2007-2008,  and  doc¬ 
umenting  the  3D  structure  of  the  heliosphere.  Consideration  and  discussion  of  these  points 
can  be  found  in  Stone  et  al.  (2005),  Heber  and  Potgieter  (2006),  Balogh  and  Erdos  (2013), 
Heber  (2013),  and  Mewaldt  (2013). 


2  Solar  Cycles  in  the  Heliosphere 

The  very  existence  and  dynamical  balance  of  the  heliosphere  are  ultimately  governed  by  the 
solar  surface  activity  which  is  driven  by  the  1 1-year  cyclic  solar  dynamo  process  leading  to 
variations  in  the  solar  wind  dynamic  pressure.  As  a  result,  the  heliospheric  extent  (distance 
to  the  termination  shock  in  the  ecliptic  plane)  would  be  maximal  on  the  decline  of  the  cycle 
when  low-latitude  high-speed  streams  are  most  prominent  (e.g.,  Richardson  and  Schwadron 
2008). 

2. 1  Solar  Wind  and  Heliospheric  Magnetic  Eield  in  the  Ecliptic  Plane 

The  temporal  variability  of  solar  magnetic  activity  (quantified  by  the  sunspot  number — 
Usoskin  2013;  Clette  et  al.  2014)  is  shown  in  Pig.  1  along  with  the  key  heliospheric 
parameters — solar  wind  speed,  proton  density,  HMP — since  1965. 

The  heliospheric  parameters  presented  here  were  measured  in-situ  at  the  Earth’s  orbit, 
i.e.  in  the  ecliptic  plane  at  1  AU.  SSN  variability  (Pig.  lA)  is  dominated  by  the  11 -year 
cycle,  varying  from  nearly  no  spots  at  solar  minimum  to  high  values  of  about  200  at  solar 
maximum.  To  date,  current  cycle  No.  24  which  started  in  December  2008  is  only  about  half 
as  large  as  the  average  of  cycles  18-23,  indicating  the  end  of  the  Modern  Grand  Maximum  of 
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1965  1970  1975  1980  1985  1990  1995  2000  2005  2010  2015 

Years 

Fig.  1  Variability  of  the  solar  and  heliospheric  parameters  for  the  period  1965-2013.  (A)  Monthly  inter¬ 
national  sunspot  numbers  (http://sidc.oma.be/).  (B)  Solar  wind  velocity  at  the  Earth’s  orbit  (OMNIweb  data 
http://omniweb.gsfc.nasa.gov/).  (C)  Solar  wind  proton  density  at  the  Earth’s  orbit  (OMNIweb  data).  (D)  He¬ 
liospheric  magnetic  field  HMF  at  the  Earth’s  orbit  (OMNIweb  data).  Red  curves  depict  the  running  13-month 
means 


solar  activity  (Solanki  et  al.  2004;  Usoskin  et  al.  2007;  Clette  et  al.  2014).  On  the  other  hand, 
heliospheric  parameters  near  Earth  show  little  variability  over  the  solar  cycle.  Solar  wind 
speed  (panel  B)  does  not  depict  a  clear  11 -year  cyclic  variability.  Richardson  et  al.  (2002) 
found  that  “For  minimum  periods,  the  Earth  is  embedded  in  high-speed  streams  ^55  %  of 
the  time  (most  prominently  on  the  decline  of  the  cycle)  versus  ^35  %  for  slow  solar  wind 
and  10  %  for  coronal  mass  ejection  (CME)-associated-structures,  while  at  solar  maximum, 
typical  percentages  are  as  follows:  high-speed  streams  ~35  %,  slow  solar  wind  ~30  %,  and 
CME-associated  ^35  Thus  the  strongest  sustained  high-speed  wind  characteristically 
occurs  on  the  decay  of  the  sunspot  curve.  Proton  density  (panel  C)  also  does  not  show  a 
signature  of  the  11 -year  cycle,  although  a  tendency  of  reduced  density  during  the  current 
weak  cycle  can  be  observed.  The  HMF  does  show  a  1 1-year  cycle  of  a  factor  of  two  as  the 
max-min  ratio.  It  is  notable  that  HMF  is  weaker  during  the  current  cycle  reaching  the  lowest 
measured  value  of  about  3.5  nT  in  2009.  The  HMF  changes  its  dipole  polarity  near  solar 
sunspot  maximum  time  every  solar  cycle  leading  to  the  22-year  magnetic  cycle  called  the 
Hale  cycle  (Oliver  2014).  For  more  detail  see  a  review  by  Owens  and  Forsyth  (2013). 
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Fig.  2  Heliolatitudinal 
distribution  of  the  daily  solar 
wind  velocities  measured  by 
Ulysses  spacecraft  (data  from 
OMNIweb).  The  upper  {blue 
dots)  and  lower  {red  dots)  panels 
are  built  for  the  solar  minimum 
(1994-1997  and  2005-2009)  and 
solar  maximum  (1991, 

2000-2003)  conditions, 
respectively.  The  plot  is  given  in 
polar  coordinates,  where  the 

radial  distance  represents  the  To 

solar  wind  speed  in  km/s,  while  E 

the  polar  angle  corresponds  to  the  ^ 

heliolatitude.  The  horizontal  % 

central  line  is  the  heliomagnetic  ^ 

equator,  while  the  grey  lines 
denote  the  range  of  the  ecliptic 
zone  in  the  heliocentric 
coordinates 


2.2  3D  View 


Since  the  ecliptic  plane  is  close  to  the  solar  equator  (slightly  tilted  with  about  7°  inclination), 
it  is  a  special  region,  and  the  heliospheric  values  measured  in  the  ecliptic  plane  may  not  be 
fully  representative  for  the  full  3D  heliosphere.  Because  of  the  tilt  between  the  ecliptic  plane 
and  the  helio-equator,  the  Earth  scans  a  narrow  range  of  heliolatitudes  (^14°)  around  the 
heliomagnetic  equator  (bounded  by  the  grey  lines  in  the  Fig.  2)  within  one  year.  However, 
this  scan  is  too  narrow  to  give  a  flavor  of  the  full  3D  variability  of  the  solar  wind.  While  HMF 
in  the  ecliptic  plane  is  roughly  representative  for  the  entire  heliosphere  (Fockwood  2013), 
the  solar  wind  velocity  varies  dramatically  at  different  helio-latitudes.  Figure  2  summarizes 
the  observations  made  onboard  the  Ulysses  spacecraft  during  its  latitudinal  scans  of  the 
heliosphere  (Wenzel  et  al.  1992).  The  upper  panel  (blue  dots)  shows  the  solar  wind  velocity 
for  the  solar  minimum  condition.  A  clear  pattern  is  observed — slow  solar  wind  dominates 
at  low  heliolatitudes,  while  fast  solar  wind  dominates  at  high  latitudes  (McComas  et  al. 
2008).  This  structure  is  caused  by  the  presence  of  large  unipolar  coronal  holes  in  the  Sun’s 
polar  regions  that  emit  fast  solar  wind.  The  situation  is  dramatically  changed  around  solar 
maxima  (lower  panel  of  Fig.  2),  when  a  mixture  of  fast  and  solar  wind  streams  appear, 
caused  by  widespread  small  coronal  holes  and  CME  activity.  Therefore,  the  solar  wind  as 
measured  in  the  ecliptic  plane  is  not  representative  for  the  full  heliosphere.  Moreover,  the 
cyclic  variability  of  the  solar  wind  velocity  at  the  ecliptic  plane  (slow  solar  wind  at  solar 
minima  and  more  frequent  fast  solar  wind  streams  at  late  declining  phase  of  the  solar  cycle) 
is  opposite  to  the  higher  latitude  pattern. 
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Fig.  3  Time  variability  of  the 
open  magnetic  flux  since  the 
Maunder  minimum  of  solar 
activity,  reconstructed  (black 
line — Vieira  and  Solanki  2010) 
and  computed  from  in  situ 
measurements  (red 
line — Lockwood  et  al.  2009) 


2.3  Open  Solar  Magnetic  Flux 

An  important  factor  for  the  cosmic  rays  in  the  heliosphere  is  the  open  magnetic  flux  (OMF) 
of  the  Sun  (Cane  et  al.  1999;  Cliver  and  Ling  2001b;  Lockwood  2013),  which  can  be  cal¬ 
culated  from  the  parameters  measured  in  situ  near  Earth  or  from  the  solar  surface  magnetic 
maps,  and  then  recalculated  to  the  OMF  at  the  source  surface  at  2.5  solar  radii  around  the 
Sun.  OMF  shows  a  great  variability  both  in  solar  cycle  and  on  the  centennial  scale  (see 
Fig.  3).  It  is  important  to  note  that  the  centennial  variability  is  great  (Lockwood  et  al.  1999; 
Solanki  et  al.  2000)  comparable  with  or  even  greater  than  the  11 -year  cycle  range  of  vari¬ 
ability.  It  is  noteworthy  that  the  recent  solar  minimum  of  2009  yielded  low  OMF  values, 
comparable  to  those  around  1900  and  during  the  Dalton  minimum  ca.  1820,  indicating  the 
end  of  the  Modern  Grand  Maximum  of  solar  activity. 

2.4  Heliospheric  Current  Sheet 

A  specific  feature  of  the  heliospheric  structure,  which  is  very  important  for  cosmic  ray  mod¬ 
ulation  (Potgieter  2013),  is  the  heliospheric  current  sheet  (HCS)  which  is  a  thin  shear  in¬ 
terface  separating  the  regions  of  the  opposite  polarity  of  HMF  emerging  from  the  Sun.  It 
corresponds  to  the  heliomagnetic  equator.  Since  the  magnetic  dipole  axis  of  the  Sun  is  tilted 
with  respect  to  rotational  axis,  this  together  with  the  Sun’s  rotation  and  radially  expanding 
solar  wind,  leads  to  formation  of  a  complicated  3D-structure,  resembling  a  ballerina’s  skirt 
(Mursula  and  Hiltula  2003;  Owens  and  Forsyth  2013).  The  waviness  of  HCS,  which  plays 
an  important  role  in  GCR  modulation  (Alanko-Huotari  et  al.  2007b),  is  defined  by  the  tilt 
angle,  which  varies  cyclically  over  the  solar  cycle  so  that  the  HCS  is  nearly  flat  around  the 
solar  minimum  and  reaches  its  maximum  values  of  above  70°  around  solar  maximum  times 
(Fig.  4).  Direct  computations  of  the  tilt  angle  started  in  mid-1970s  but  indirect  estimates 
(Pishkalo  2006)  can  be  done  even  before  that  based  on  observations  of  solar  eclipses  since 
1870  (Fig.  4).  The  results  suggest  that  the  cyclic  variability  of  the  HCS  warp  does  not  de¬ 
pend  on  the  strength  of  the  solar  cycle — the  tilt  angle  varies  to  its  full  extent  for  both  high  (in 
the  second  half  of  20th  century)  and  low  (ca.  1900)  solar  cycles.  It  does  not  depict  a  centen¬ 
nial  variation.  Accordingly,  we  may  expect  a  similar  pattern,  for  example  for  the  Maunder 
minimum  condition. 

2.5  Grand  Minima  of  Solar  Activity 

Of  special  interest  is  a  question  what  are  the  heliospheric  conditions  during  a  grand  mini¬ 
mum  of  solar  activity,  the  most  famous  example  being  the  Maunder  minimum  in  1645-1700 
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Fig.  4  Variability  of  the  HCS  tilt 
angle  as  a  function  of  the  1 1-year 
solar  cycle  phase.  The  red  dots 
represent  the  HCS  tilt  angle 
(radial  model)  averaged  over 
cycles  21-23  (Wilcox  Solar 
Observatory  http://wso. Stanford. 
edu/Tilts.html),  the  blue  stars 
depict  the  tilt  angle  reconstructed 
from  solar  eclipses  since  1870 
(Pishkalo  2006),  and  the  blue 
curve  is  a  7-point  running  mean 
over  the  blue  stars  (cf. 
Alanko-Huotari  et  al.  2007a) 


Solar  cycle  phase 


with  almost  no  sunspots  on  the  solar  surface  (Eddy  1976).  Despite  the  lack  of  sunspots  a 
weak  cyclic  variability  was  observed  during  the  Maunder  Minimum  in  geomagnetic  and 
cosmic  ray  indices  (Beer  et  al.  1998;  Usoskin  et  al.  2001;  Miyahara  et  al.  2004;  Berggren 
et  al.  2009;  Owens  et  al.  2012)  suggesting  that  the  solar  magnetic  cycle  was  still  operating 
in  the  heliosphere.  Different  estimates  of  the  HMF  and  solar  wind  parameters  exist,  but  the 
current  paradigm  (Oliver  et  al.  1998;  McCracken  2007a;  Steinhilber  et  al.  2010;  Owens  and 
Forsyth  2013;  Fockwood  and  Owens  2014)  is  that  the  HMF  was  weak,  about  2  nT,  and  solar 
wind  in  the  ecliptic  plane  remained  slow  as  emanating  from  hot  parts  of  the  corona  even  in 
the  absence  of  active  sunspot  regions.  The  likelihood  that  the  solar  wind  was  slow  during  the 
Maunder  Minimum  was  underscored  by  the  recent  solar  cycle  minimum  between  cycle  23 
and  24,  when  Earth  was  embedded  in  slow  solar  wind  ~70  %  of  the  time  in  2009  (Oliver  and 
Fing  2011).  However,  because  of  the  imbalance  between  the  different  processes,  the  phase 
relation  between  OMF  (and  hence  cosmic  rays)  and  the  solar  magnetic  cycle  was  probably 
inverted  with  respect  to  normal  solar  cycles  (Owens  et  al.  2012).  Such  grand  minima  occur 
every  now  and  then  (Usoskin  2013)  and  correspond  to  a  special  mode  of  the  solar  dynamo 
(Usoskin  et  al.  2014). 


3  Solar  Cycles  in  Galactic  Cosmic  Rays  as  Observed  at  Earth 

In  this  section  we  discuss  how  the  conditions  in  the  heliosphere  are  reflected  by  the  results 
of  cosmic  ray  observations.  Galactic  cosmic  rays  are  highly  energetic  fully  ionized  nuclei 
coming  into  the  heliosphere  from  outside  the  solar  system  with  a  roughly  constant  flux. 
GCRs  with  kinetic  energy  in  the  range  from  a  few  hundred  MeV  up  to  100  GeV  can  be 
regarded  as  heliospheric  probes.  Transport  of  cosmic  rays  in  the  heliosphere  is  governed  by 
several  processes,  viz.  diffusion  due  to  scattering  on  the  HMF  irregularities,  convection  with 
the  radially  expanding  solar  wind  plasma,  drift  in  the  large-scale  heliospheric  magnetic  field 
including  that  along  the  HCS,  and  adiabatic  energy  losses  due  to  the  divergence  of  solar 
wind  (see  Potgieter  2013,  for  details).  All  these  heliospheric  characteristics  depend  on  the 
phase  of  the  solar  magnetic  cycle  and  on  the  solar  activity  level.  As  a  result,  the  GCR  flux 
is  modulated  by  the  solar  magnetic  activity.  Therefore,  by  measuring  GCR  flux  at  Earth  one 
can  study  the  1 1-year  and  22-year  cycles  in  the  heliosphere. 

Figure  5  presents  time  profiles  of  SSN  as  well  as  GCR  intensity  as  observed  by  different 
instruments  sensitive  to  different  energy  ranges  of  primary  particles.  The  GCR  records  vary 
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Years 

Fig.  5  Time  variability  of:  (A)  sunspot  numbers  (http://sidc.oma.be/sunspot-data/);  (B)  cosmic  ray  relative 
intensity  measured  in  space  onboard  IMPS  and  PAMELA  spacecraft;  (C)  cosmic  ray  relative  intensity  mea¬ 
sured  by  high  latitude  balloons  (LPI);  (D)  cosmic  ray  relative  intensity  measured  by  ground-based  neutron 
monitors,  polar  Oulu  and  equatorial  Huancayo/Haleakala.  Relative  variations  of  cosmic  rays  are  normalized 
to  100  %  in  March  1987 


cyclically  in  rough  anti-phase  to  SSN  with  a  delay  which  depends  on  the  direction  of  the 
CR  drifts  in  the  heliosphere  (Usoskin  et  al.  1998;  Cliver  and  Ling  2001a).  The  depth  of  the 
modulation  depends  greatly  on  the  CR  energy;  it  varies  from  a  few  %  for  an  equatorial  NM 
to  ^20  %  for  a  polar  NM  (see  panel  D);  a  factor  of  two  for  stratospheric  balloon-borne  data 
(panel  C);  and  up  to  an  order  of  magnitude  for  space-borne  data  (panel  B).  One  can  also  see 
alternation  of  sharp  and  flat  peaks  in  the  GCR  intensity  that  is  related  to  the  dominant  drifts 
in  the  heliosphere  (Potgieter  2013).  An  important  feature  is  that  the  GCR  intensity  in  2008- 
2010  was  the  record  high  for  the  entire  series  reflecting  the  weak  solar  activity.  The  GCR 
intensities  (Fig.  5)  follow  the  (inverted)  1 1 -year  solar  cycle  much  better  than  any  individual 
heliospheric  plasma  and  magnetic  field  parameter  observed  in  the  near  Earth  space  (Fig.  1). 
On  the  other  hand,  GCR  flux  is  well  related  to  the  OMF  from  coronal  holes  at  latitudes 
<45°  (Cliver  and  Ling  2001b;  Alanko-Huotari  et  al.  2006).  This  implies  that  the  local  in 
situ  heliospheric  measurements  at  1  AU  in  the  ecliptic  plane  may  not  necessarily  represent 
the  global  heliospheric  changes,  since  the  GCR  modulation  is  affected  by  the  plasma  and 
magnetic  field  parameters  in  the  whole  heliosphere  rather  than  locally. 
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3.1  Cosmic  Ray  Detectors  and  Their  Energy  Responses 

In  order  to  study  the  1 1 -year  cycle  in  CR  intensity  one  needs  long-term  time  series  of  mea¬ 
surements  of  the  charged  particle  flux.  Since  the  GCR  modulation  depends  on  particle  en¬ 
ergy,  detectors  sensitive  to  different  energy  ranges  of  CR  are  needed.  The  low  energy  part  of 
the  GCR  spectrum,  up  to  tens  of  GeV,  is  affected  by  the  heliospheric  modulation,  but  there 
is  an  indication  that  a  weak  modulation  may  affect  even  ^100  GeV  particles  (Potgieter  and 
Strauss  2013). 

In  the  energy  range  of  0.5-15  GeV,  the  geomagnetic  field  and  the  Earth’s  atmosphere 
serve  as  an  energy  separator  (see  e.g.,  Grieder  2001).  Eirst  attempts  to  initiate  a  permanent 
ground-based  regular  CR  observational  network,  with  identical  ionization  chambers  (IC), 
refer  back  to  1930s  (Eorbush  1954).  ICs  at  sea  level  detect  secondary  muons  generated  in 
the  atmosphere  by  protons  with  energies  >4  GeV  impinging  at  the  top  of  the  atmosphere. 

However,  a  real  impetus  to  regular  CR  monitoring  was  given  by  John  Simpson  with 
invention  of  a  neutron  monitor  (Simpson  1958,  2000).  A  world-wide  network  of  ground- 
based  neutron  monitors  (NMs)  was  developed  in  1950s  and  included,  in  different  times, 
more  than  a  hundred  NMs,  with  more  than  40  NMs  operating  presently  (Moraal  et  al.  2000; 
Mavromichalaki  et  al.  2011).  The  NM  data  set  is  available  from  mid-1950s  up  to  now  (see, 
e.g.,  Eig.  5d).  Because  of  the  geomagnetic  shielding,  each  NM  accepts  particles  above  some 
cut-off  energy /rigidity  defined  by  the  location  of  the  NM  and  the  geomagnetic  field  strength. 
The  shielding  is  absent  in  polar  and  maximum  in  equatorial  regions  where  the  cut-off  rigidity 
is  about  15  GV.  Thus,  the  global  NM  network  can  be  considered  as  a  single  device  allowing 
to  obtain  information  on  the  CR  particle  fluxes,  their  energy  and  anisotropy.  Reconstructions 
of  the  GCR  flux  based  on  the  NM  network  data  agree  well  with  the  direct  balloon-  and  space- 
borne  measurements,  making  the  NM  network  a  very  useful  instrument  to  monitor  CR  flux 
(Usoskin  et  al.  2011). 

The  problem  of  converting  the  data  obtained  on  the  ground  and  in  the  atmosphere 
into  characteristics  of  the  primary  CRs  requires  knowledge  of  the  response/yield  function 
of  different  devices  and  understanding  of  the  atmospheric/magnetospheric  transport  (e.g., 
Caballero-Lopez  and  Moraal  2012).  As  a  measure  of  the  sensitivity  of  a  NM  to  CR  energy, 
Alanko  et  al.  (2003)  introduced  a  concept  of  the  effective  energy  for  a  NM,  £eff,  such  that 
the  count  rate  of  a  given  NM  is  proportional  to  the  flux  of  GCR  with  energy  above  E^eff, 
which  varies  between  6.5  GeV  for  a  polar  NM  to  30-40  GeV  for  an  equatorial  NM. 

Another  common  type  of  ground-based  cosmic-ray  detector  is  a  muon  telescope  mea¬ 
suring  the  muon  component  of  the  cosmic-ray  induced  atmospheric  cascade.  They  typically 
have  higher  effective/median  energy  (40-60  GeV — see  Jamsen  et  al.  2007)  than  NMs. 

The  ground-based  monitoring  is  complemented  by  CR  observations  in  the  Earth’s  at¬ 
mosphere.  During  1951-1969,  Neher  (1967,  1971)  accomplished  a  series  of  regular  but 
infrequent  balloon  measurements  of  the  cosmic  ray  induced  ionization.  In  mid-1950s,  a 
Soviet/Russian  program  (Charakhchyan  1964;  Stozhkov  et  al.  2009)  lead  by  the  Lebedev 
Physical  Institute  (LPI)  of  long-term  balloon-borne  observations  of  the  CR  fluxes  at  several 
latitudes  in  the  Earth’s  atmosphere  had  started  and  continues  till  now  (see  Eig.  5c).  Balloons 
were  launched  on  the  daily  basis  until  1990,  and  three  times  a  week  since  then.  The  LPI 
balloon  measurements  get  data  from  the  ground  level  up  to  the  altitude  of  about  30  km.  This 
data  corresponds  to  £eff  ^  1-2  GeV  (Bazilevskaya  et  al.  2013). 

Numerous  balloon-borne  detectors  have  been  launched  to  measure  the  CR  spectrum  at 
shallow  atmospheric  depth,  such  as,  e.g.,  BESS  instrument  (e.g.  Shikaze  et  al.  2007).  How¬ 
ever,  such  measurements,  while  providing  direct  measurements  of  the  GCR  spectrum,  are 
campaign-based  and  do  not  provide  long-term  monitoring. 
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Cosmic  rays  with  energy  below  several  hundred  MeV  can  be  recorded  only  by  spacecraft- 
borne  detectors  because  such  particles  are  absorbed  in  the  upper  atmosphere.  Only  a  few  lim¬ 
ited  data  sets  of  sufficient  duration  are  available,  among  them  the  results  from  the  IMPS,  and 
the  ACE  spacecraft  (McDonald  1998;  McDonald  et  al.  2010),  and  the  more  recent  PAMELA 
(Adrian!  et  al.  2013)  and  AMS  (Ting  2013)  orbital  spectrometers.  A  record  is  shown  in 
Eig.  5b  and  corresponds  to  the  flux  of  protons  with  energy  120-230  MeV  (Bazilevskaya 
et  al.  2013). 

3.2  CR  Transport  in  the  Heliosphere 

Temporal  variability  of  GCR  intensity  in  the  energy  range  affected  by  solar  modulation  pro¬ 
vides  indirect  information  on  the  state  of  the  whole  heliosphere,  not  only  in  the  vicinity 
of  Earth.  Eor  this  purpose,  one  needs  to  understand  the  processes  driving  the  CR  trans¬ 
port  in  the  heliosphere.  The  basics  of  the  CR  transport  equation  were  introduced  in  the 
1960s  (Parker  1965;  Jokipii  1966;  Dolginov  and  Toptygin  1967;  Gleeson  and  Axford  1967; 
Krymskij  1969)  and  further  developed  by  a  number  of  groups  (e.g.,  Dorman  2006).  Mod¬ 
ern  CR  transport  models  are  sophisticated  and  precise  (for  a  review  see  Potgieter  2013,  and 
references  therein),  allowing  to  study  numerically  different  effects  and  processes.  In  order 
to  study  the  full  process,  however,  one  needs  detailed  information  on  the  local  interstellar 
spectrum  of  GCRs,  as  well  as  geometrical  structure,  polarity,  strength,  and  turbulence  level 
of  the  HME  and  the  solar- wind  speed  throughout  of  the  3D  heliosphere.  Unfortunately,  as 
discussed  above,  there  are,  with  rare  exceptions,  only  local  near-Earth  observations  of  the 
HME  parameters.  Therefore,  the  models  have  to  use  assumptions,  e.g.,  regarding  elements 
of  the  diffusion  tensor  that  cannot  be  measured  directly.  On  the  other  hand,  modulation  of 
GCRs  with  energy  above  several  hundred  MeV  can  be  formally  described  by  the  so-called 
force-field  model  (e.g.,  Gleeson  and  Axford  1968;  Caballero-Lopez  and  Moraal  2004).  This 
model  parameterizes  the  GCR  spectrum  with  a  single  formal  parameter,  the  modulation  po¬ 
tential.  The  momentary  value  of  the  modulation  potential  can  be  calculated  from  the  data 
of  the  ground-based  NM  network  (Usoskin  et  al.  2011).  On  the  long-term  timescale,  the 
modulation  parameter  can  be  reconstructed  from  data  of  cosmogenic  isotopes  ^"^C  and  ^^Be 
measured  in  terrestrial  archives  (see,  e.g..  Beer  et  al.  2012;  Usoskin  2013). 

It  has  been  demonstrated  that  this  approach,  while  not  pretending  to  be  a  physical  in¬ 
terpretation,  provides  a  very  useful  single-parameter  parameterization  of  the  GCR  energy 
spectrum  near  Earth  (see  formalism  in  Usoskin  et  al.  2005),  especially  for  long-term  stud¬ 
ies.  Nonetheless,  a  more  sophisticated  approach  is  needed  to  give  insight  into  detail  of  CR 
interaction  with  the  heliospheric  magnetic  field. 

3.3  CR  Diagnostic  of  the  Heliosphere 

Modulation  of  GCR  in  the  heliosphere  is  largely  defined  by  the  turbulent  component  of  HME 
responsible  for  particle  diffusion,  as  the  mean  free  path  depends  on  the  power  spectrum  of 
the  magnetic  irregularities  in  the  heliosphere  (Jokipii  1966).  Recent  models  (e.g.,  Eerreira 
and  Potgieter  2004)  are  good  in  explaining  the  temporal  variability  of  the  observed  GCR 
fluxes,  implying  that  the  basic  conception  of  CR  modulation  is  correct.  However,  the  value 
of  the  mean  free  path  remains  unknown  because  of  the  complicated  structure  of  the  turbulent 
interplanetary  magnetic  field  (Giacalone  2013).  The  HME  has,  in  the  rough  approximation, 
a  form  of  an  Archimedian  spiral  (Parker  1965)  that,  leads  to  gradient  and  curvature  drifts  of 
CR  in  the  interplanetary  medium.  Account  for  motion  of  the  magnetic  field  foot-points  on 
the  Sun  leads  to  a  more  realistic  Eisk-type  HME  (Eisk  1996)  which  is  particularly  important 
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at  high  latitudes.  Important  is  also  a  drift  at  the  HCS  which  is  an  extreme  case  of  the  gradi¬ 
ent  drift.  The  role  of  the  drift  effects  in  CR  modulation  was  pointed  out  by  Jokipii  and  Levy 
(1977)  and  intensively  studied  since  then  (e.g.,  Alanko-Huotari  et  al.  2007b;  Strauss  et  al. 
2012b).  When  the  solar  magnetic  field  is  directed  outward  from  the  Sun  in  the  north  polar 
region  (the  so-called  A  >  0  periods),  positively-charged  particles  drift  inwards  in  the  helio¬ 
spheric  polar  regions  and  drift  outwards  along  the  HCS.  This  facilitates  CR  access  into  the 
inner  heliosphere  leading  to  the  long  fiat  shape  of  CR  intensity  maxima  as,  e.g.,  ca.  1976  and 
1996.  The  drift  has  the  opposite  pattern  during  A  <  0  periods,  forming  favorable  conditions 
for  the  CR  access  to  the  inner  heliosphere  when  the  HCS  is  fiat.  This  leads  to  sharp  high 
peaks  of  CR  intensity  in  1965,  1987  and  2009.  This  drift-dominated  modulation  is  typical 
(but  see  the  next  section)  for  minima  of  solar  activity  cycles  when  the  heliospheric  structure 
is  fairly  regular.  The  shape  alternation  of  CR  maxima  is  seen  in  Fig.  5  as  a  prominent  signa¬ 
ture  of  the  drift  effect  in  the  CR  modulation.  However,  during  periods  of  high  solar  activity 
drifts  are  less  important  (Krainev  and  Kalinin  2013)  since  the  modulation  is  mostly  driven 
by  diffusion  and  propagating  barriers  (Burlaga  et  al.  1985). 

3.4  Unusual  Minimum  of  Solar  Activity  Between  Solar  Cycles  23  and  24 

As  one  can  see  in  Fig.  5,  the  intensity  of  GCR  reached  in  2009  its  record  highest  value,  over 
the  entire  period  of  direct  observations  since  1950s  (see,  e.g.,  Heber  et  al.  2009;  Ahluwalia 
et  al.  2010;  McDonald  et  al.  2010;  Moraal  and  Stoker  2010;  Krainev  et  al.  2013).  This  hap¬ 
pened  during  the  minimum  of  solar  activity  between  solar  cycles  23  and  24,  which  is  known 
to  be  the  weakest  heliospheric  activity  period  for  the  space  era  (Gibson  et  al.  2011;  McCo- 
mas  et  al.  2013).  This  period  is  characterized  by  the  record-low  HMF  intensities,  reduced 
HMF  turbulence,  change  in  the  inhomogeneity  spectrum  of  the  HMF,  and  reduced  solar- 
wind  dynamic  pressure.  The  GCR  intensity  reached  a  sharp  maximum  in  2009  following 
the  flattening  of  the  HCS  and  then  quickly  dropped  to  the  normal  level  after  the  sudden 
increase  of  the  HCS  tilt  angle  in  early  2010  (McComas  et  al.  2013).  The  GCR  intensity 
enhancement  was  observed  in  a  wide  range  of  particle  energies  (Fig.  5).  It  was  as  great  as 
60  %  at  a  few  hundred  MeV  energy  (panel  B)  and  only  a  few  %  in  the  data  of  a  polar  NM 
(panel  D).  On  the  other  hand,  the  CR  intensity  as  recorded  at  mid-  and  low-latitudes  (geo¬ 
magnetic  cutoff  rigidity  above  4  GV)  did  not  show  any  excess  in  2009  compared  to  1987 
(Moraal  and  Stoker  2010).  The  unusual  rigidity  dependence  of  the  GCR  increase  suggests 
that  significant  changes  took  place  in  the  heliosphere  properties  during  the  recent  solar  cy¬ 
cle  minimum  (Mewaldt  et  al.  2010;  Bazilevskaya  et  al.  2012;  Krainev  and  Kalinin  2013; 
Potgieter  et  al.  2014).  These  facts  indicate  that  mainly  lower  energy  GCRs  (<10  GeV)  were 
affected  by  the  favorable  heliospheric  conditions  but  not  higher  energies  (see,  e.g.  Usoskin 
et  al.  201 1).  Cliver  et  al.  (2013)  pointed  out  that  the  record  increase  in  cosmic  ray  intensity 
in  2009 — which  was  accompanied  by  a  decrease  in  HMF  and  an  increase  in  the  tilt  angle 
relative  to  the  previous  minimum — challenged  the  dominance  of  drift  effects  for  modulation 
at  1 1-yr  minima  of  the  solar  cycle.  More  recently,  Potgieter  et  al.  (2014)  concluded  that  “the 
2009  modulation  minimum  could  be  described  as  more  ‘diffusion  dominated’  than  previous 
solar  minima”  and  that  “diffusion  contributed  ~50  %  of  the  total  cosmic  proton  intensities 
observed  at  Earth  while  particle  drifts  contributed  the  other  50  This  implies  that  the 
scaling  of  the  HMF  turbulence  changed  during  the  recent  solar  cycle  minimum  compared  to 
the  previous  minimum  period,  although  a  more  detailed  analysis  and  modelling  is  required 
to  fully  understand  this  process  (Strauss  et  al.  2012a). 
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4  Cycles  in  Solar  Energetic  Particles 

4. 1  Solar  Proton  Events 

Unlike  many  solar  and  geophysical  measurements  such  as  SSN  or  geomagnetic  indices,  the 
history  of  identification  of  solar  proton  events  (SPEs)  is  not  as  long.  The  first  SPEs  were 
directly  recorded  by  ICs  in  Cheltenham  (USA),  Godhavn  (Greenland),  and  Christchurch 
(New  Zealand)  in  1942  as  clear  increases  on  28-Eeb  and  7-Mar  (Lange  and  Eorbush  1942). 
However,  these  increases  were  not  ascribed  to  particles  from  the  Sun  until  the  third  increase 
was  observed  at  these  same  locations  on  25-Jul-1946  (Eorbush  1946).  A  fourth  increase 
(19-NOV-1949)  was  also  recorded  by  other  detectors,  e.g.,  by  an  IC  in  Climax,  USA  (Eor¬ 
bush  et  al.  1950),  muon  detectors  at  several  locations  and  a  prototype  NM  at  Manchester, 
UK  (Adams  1950).  While  direct  measurements  of  SPEs  prior  to  the  19th  solar  cycle  were 
limited  to  the  ICs  and  early  NMs,  there  were  indirect  measurements  from  vertical  incidence 
ionospheric  soundings  made  at  high  latitude  locations  starting  in  1938.  In  a  comparison 
of  polar  cap  absorption  events  in  1956-1959  with  records  of  vertical  incidence  soundings, 
Svestka  (1966)  compiled  a  fairly  homogeneous  set  of  data  on  strong  polar  cap  absorption 
events  for  1938-1959.  During  the  19th  solar  cycle,  lower  energy  SPEs  could  be  identified 
by  ionospheric  measurements,  and  instruments  on  balloons  and  early  rocket  and  satellite 
experiments,  that  were  able  to  detect  solar  protons  at  lower  energies  (MeV  range).  The  solar 
proton  catalog  edited  by  Svestka  and  Simon  (1975)  provides  a  compendium  of  SPEs  of  the 
19th  solar  cycle,  which  is  not,  however,  intended  to  provide  a  homogeneous  database;  nev¬ 
ertheless  it  is  the  best  data  available  during  that  time  period.  Similar  catalogs  were  compiled 
also  later  under  the  supervision  of  Yu.I.  Logachev^  spanning  the  period  of  for  1970-1996. 

The  installation  of  the  standardized  NMs  starting  in  1951  considerably  enhanced  the 
ability  to  detect  large  relativistic  SPEs.  The  event  of  23-Eeb-1956  was  the  largest  increase 
recorded  by  ground-based  NMs  and  muon  detectors;  the  Leeds  NM  recorded  an  increase  of 
4554  %  over  a  15-minute  interval,  and  the  Moscow  muon  detector  recorded  an  increase  of 
300  %.  Throughout  the  19th  solar  cycle  (1954-1965)  ground-based  NMs  and  muon  detec¬ 
tors  were  essentially  the  only  instruments  that  could  identify  high  energy  relativistic  SPEs. 
Accordingly,  these  events  were  called  ground-level  enhancements,  or  ground-level  events 
(GLEs). 

Starting  with  the  IMP  satellite  program  in  1964,  it  became  possible  to  assemble  an  al¬ 
most  homogeneous  database  of  SPEs  recorded  at  Earth.  In  1976,  the  NOAA/USAE  Space 
Environment  Services  Center  in  Boulder,  Colorado  defined  a  significant  solar  proton  event 
(SEP)  as  any  event  with  a  proton  flux  greater  than  10  protons  (cm^  s  sr)“^  above  10  MeV. 
Using  this  criteria.  Shea  and  Smart  (1990)  compiled  a  list  of  SPEs  from  1955  until  1986 
(i.e.  solar  cycles  19-21).  While  an  update  of  that  list  is  beyond  the  scope  of  this  paper,  the 
events  in  solar  cycles  22  through  cycle  24  (first  5  years)  have  been  included  in  the  following 
statistics  and  graphs. 

4.1.1  Statistics  of  Solar  Proton  Events 

Using  the  NOAA  criteria  of  identifying  a  SPE  (see  above)  we  have  identified  all  events 
that  meet  this  criteria  in  solar  cycles  19  through  year  5  in  solar  cycle  24.  Since  the  list  of 
SPEs  maintained  by  NOAA  does  not  identify  independent  injections  of  solar  protons  within 


^http://www.  wdcb.ru/stp/data/PRCATFINAL/. 
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Table  1  Summary  of  Solar  Proton  Events  for  Solar  Cycles  19  through  the  first  five  months  of  Solar  Cycle  24, 
including:  start  months  and  duration  (in  months)  of  solar  cycles,  the  total  number  of  discrete  events,  iVsPE^ 
total  number  of  GLEs  iVoLE^  percentage  of  GEE  Pgle  vs.  the  SPEs,  and  the  omnidirectional  fluence  (in 
10^^  cm“^)  of  >10  MeV  protons,  Fio 


Cycle 

Start 

Duration 

VsPE 

Vgle 

Pole 

Go 

19 

May  1954 

126 

65 

10 

15.4 

7.2 

20 

Nov  1964 

140 

72 

13 

18.0 

2.2 

21 

Jul 1976 

123 

81 

12 

14.8 

1.8 

22 

Oct  1986 

120 

84 

15 

17.8 

5.8 

23 

Oct  1996 

146 

103 

16 

15.5 

7.5 

24 

Dec  2008 

60 

29 

1 

3.4 

1.2 

a  major  episode  of  activity  unless  the  >10  MeV  flux  falls  below  the  event  threshold,  we 
have  examined  each  of  the  events  for  independent  and  additional  sources  of  solar  protons. 
Increases  associated  with  the  arrival  of  interplanetary  shocks  at  Earth  were  not  assumed 
to  be  separate  events  but  considered  part  of  the  previous  identified  event.  Occasionally  a 
proton-producing  region  of  activity  on  the  far  Eastern  limb  of  the  Sun  can  be  associated 
with  a  small  proton  enhancement  that  does  not  reach  the  NOAA  proton  event  criteria.  When 
the  interplanetary  shock  from  the  coronal  mass  ejection  (CME)  of  such  an  event  arrives 
at  Earth,  the  additional  protons  accelerated  by  the  shock  may  add  to  the  already  enhanced 
flux  such  that  the  SPE  criteria  will  be  met.  While  we  considered  this  as  a  shock  associated 
increase,  it  will  be  listed  as  a  proton  event  on  the  NOAA  listing.  An  example  of  this  type  of 
event  occurred  on  23-Jun-2013. 

The  events  we  have  identified  as  separate  and  unique  solar  proton  events  were  tallied  in 
12-month  arrays  starting  with  the  first  month  of  each  solar  cycle.  The  start  of  each  solar 
cycle  was  taken  as  the  month  after  the  minimum  smoothed  sunspot  number  (see  Table  1). 

The  total  omnidirectional  fluence  (time-integrated  flux)  of  >  10  MeV  protons  near  Earth, 
denoted  as  Fio,  is  often  considered  as  a  characteristic  integral  parameter  of  SPEs.  The  values 
of  Fio  for  the  past  6  solar  cycles  (cycle  24  incomplete  yet)  are  given  in  Table  1.  These 
values  differ  by  factors  of  4  from  cycle  to  cycle,  which  can  be  attributed  to  the  location  of 
the  parent  solar  activity  associated  with  the  various  increases  over  a  solar  cycle  (Shea  and 
Smart  2012).  Eor  example  one  strong  event,  e.g.,  in  August  1972  can  completely  dominate 
the  total  fluence  for  the  entire  solar  cycle  (Shea  and  Smart  1990). 

Eigure  6  illustrates  the  distribution  of  SPEs  throughout  solar  cycles  19-23  and  the  first 
five  years  of  cycle  24.  The  distribution  of  events  identified  by  Svestka  (1966)  for  cycles  17- 
18  are  shown  in  Eig.  7  using  the  same  convention.  While  caution  should  be  used  in  utilizing 
these  data  for  comparison  with  the  modern  SPE  identification,  the  distributions  shown  in 
Eig.  7  are  similar  to  those  shown  in  Eig.  6.  Eigure  8  presents  the  average  number  of  SPEs 
over  solar  cycles  19-23  together  with  the  average  sunspot  number  summed  over  the  same 
time  period.  While  no  coherent  pattern  of  the  SPE  occurrence  is  seen  in  Eig.  6  for  individual 
cycles,  the  overall  pattern  of  the  SPE  number  appears  to  be  consistent  with  the  average 
sunspot  number  when  summed  over  five  solar  cycles.  The  distribution  of  GLEs  over  solar 
cycles  19-24  is  shown  in  Eig.  9.  In  comparing  it  with  Eig.  8,  it  appears  that  the  occurrence 
of  GLEs  is  more  widely  distributed  over  the  solar  cycle  than  the  >10  MeV  proton  events. 
This  may  either  reflect  a  physical  difference  or  low  statistics  of  GLEs  in  this  time  period 
(67  GLEs  vs.  405  SPEs).  When  the  four  GLEs  of  solar  cycles  17  and  18,  observed  by  muon 
detectors,  are  included,  this  distribution  becomes  even  wider. 
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Fig.  6  Distribution  of  SPEs  >10  MeV  over  each  solar  cycle  19-24  (see  Table  1  for  dates).  These  are  shown 
in  12-month  increments  starting  with  the  year  after  the  sunspot  minimum  for  each  cycle.  Only  the  first  five 
years  of  cycle  24  (i.e.  through  November  2013)  are  shown.  The  histograms  are  the  yearly  sunspot  number 
(divided  by  10)  summed  for  each  of  the  same  time  periods 


Fig.  7  The  distribution  of  SPEs  as  identified  by  Svestka  (1966)  for  solar  cycles  17  and  18.  The  number  of 
events  have  been  binned  in  12-month  intervals  from  the  start  of  cycle  17  (October  1933)  and  the  start  of  cycle 
18  (March  1944).  The  histograms  are  the  yearly  sunspot  number  (divided  by  10)  summed  for  each  of  the 
same  time  periods 


Fig.  8  Left  side:  Average  number  of  SPEs  each  year  over  five  solar  cycles  (1954-2008).  Right  side:  Average 
sunspot  number  (divided  by  10)  over  five  solar  cycles  (1954-2008).  The  yearly  values  have  been  calculated 
for  each  12-month  period  starting  with  the  first  month  after  the  minimum  value  in  the  smoothed  sunspot 
number 
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Fig.  9  Occurrence  frequency  of 
GLEs  over  a  representative 
“11-year”  solar  cycle.  Each 
complete  cycle  since  1942  has 
been  divided  into  1 1  equal 
intervals  and  the  number  of 
GLEs  has  been  distributed  over 
the  1 1  intervals  and  summed 
accordingly.  The  black  shading 
indicates  the  six  GLEs  with 
7^200  >  3.0  X  10^  proton/cm^ 
from  1956-present  (see  Sect.  4.3) 


Cycles  18-23 


1  23456789  10  11 


Solar  Cycle  Interval 


4.1.2  Ground  Level  Enhancements 

GLEs  are  relatively  rare  events  with  an  almost  homogeneous  database  since  1956,  and  as 
such,  have  been  widely  studied.  Until  2014,  there  have  been  71  GLEs  since  the  first  events 
identified  in  1942  with  a  mean  rate  of  about  one  event  per  year  (e.g.,  Cliver  et  al.  1982; 
Cliver  2006;  Gopalswamy  et  al.  2012,  2013)  but  unevenly  distributed  over  the  solar  cycles 
(Eig.  10).  Ten  GLEs  were  identified  in  the  19th  solar  cycle,  but  a  few  small  events  may  have 
been  missed  because  of  the  anisotropy  (not  known  at  that  time),  the  sparsity  of  detectors 
particularly  in  the  polar  regions,  and  a  low  flux  of  protons  >450  MeV  to  register  as  an 
increase  on  a  NM.  In  the  mid  1970s,  the  cosmic  ray  community  adopted  a  criterion  for  an 
event  to  be  identified  as  a  GLE,^  that  is  two  independent  NMs  have  recorded  a  statistically 
significant  increase.  We  note  that  an  enhancement  of  06-Jan-2014  observed  at  the  South 
Pole  NM  (Thakur  et  al.  2014)  was  a  not  a  GLE  according  to  the  formal  criterion.  As  shown 
in  Eig.  10,  these  events  can  occur  at  any  time  of  the  solar  cycle  but  tend  to  be  on  the  rising 
and  decreasing  parts  of  the  cycle.  With  the  exception  of  solar  cycle  22,  there  is  a  distinct  gap 
near  sunspot  maximum  indicative  of  the  Gnevyshev  Gap  (Gnevyshev  1967,  1977). 

One  important  aspect  of  GLEs  is  the  degree  of  anisotropy  that  may  be  present,  particu¬ 
larly  for  events  resulting  from  solar  activity  on  the  western  solar  hemisphere  (see  Sect.  3.2 
of  Shea  and  Smart  2012).  In  addition  to  the  normal  anisotropy  typically  present  during  the 
first  20-30  minutes  of  many  GLEs,  McCracken  et  al.  (2012)  have  identified  high-energy 
impulsive  GLEs  that  constitute  the  first  phase  of  some  very  large  GLEs.  These  authors  of¬ 
fered  several  constraints  that  must  be  met  by  any  putative  acceleration  mechanism  for  these 
unusual  events.  A  comprehensive  study  of  GLEs  needs  to  address  the  degree  of  anisotropy 
as  well  as  other  defining  factors  such  as  the  type  of  detector  recording  the  event,  and  the 
geomagnetic  cutoff  rigidity  and  altitude  of  the  detector  (e.g.,  Mishev  et  al.  2014). 

While  the  emphasis  on  the  detection  of  GLEs  since  1955  has  been  on  NM  measurements, 
the  importance  of  the  IC  and  muon  measurements  should  not  be  ignored  even  though  the 
particle  energies  in  most  GLEs  do  not  have  sufficient  flux  above  the  4  GeV  threshold  for 
detection  by  muon  detectors.  Of  the  67  GLEs  between  1955  and  2014,  only  5  have  been 
recorded  by  the  world  wide  network  of  muon  detectors  as  illustrated  in  Eig.  11.  This  makes 


^The  official  master  list  of  GLEs  is  available  at  http://gle.oulu.fi. 
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Fig.  10  Timing  of  the  7 1  GLEs 
observed  from  1942  to  the 
present  (grey  bars)  along  with 
sunspot  numbers.  Solar  cycles 
are  separated  by  dashed  lines 


Year 


Fig.  11  GLEs  recorded  my 
muon  detectors  (given  in 
percentage  increase)  at  the  sea 
level  from  1942  to  2010.  The 
1940  events  are  bi-hourly  data. 
The  1950  event  are  hourly  data 
and  the  events  from  1960  to  2010 
are  5  minute  data 


the  first  four  GLEs  in  solar  cycles  17  and  18  very  impressive  as  massive  events.  The  asso¬ 
ciated  solar  regions  for  the  events  of  28-Feb-1942  and  25-Jul-1946  were  located  at  ^4  °E 
and  15  °E  respectively  (Duggal  1979).  The  event  on  14-Jul-2000  was  also  associated  with 
central  meridian  activity  (7  °W)  but  the  muon  detector  increase  was  an  order  of  magnitude 
lower  than  for  the  1946  event.  The  remaining  events  recorded  by  muon  detectors  were  as¬ 
sociated  with  solar  activity  on  the  far  western  side  of  the  Sun.  It  should  be  noted  that  the 
increases  for  the  events  in  the  1940s  were  from  bi-hourly  records;  hourly  data  would  have 
more  likely  been  considerably  higher. 

4.2  Coronal  Mass  Ejections  and  SEP  Events 

It  is  well-accepted  that  the  free  energy  for  solar  activity  resides  in  solar  magnetic  fields, 
specifically  in  the  low  coronal  fields  of  active  regions  (e.g.,  Hudson  2011).  Thus,  it  is  not 
surprising  that  the  source  active  regions  of  SPEs  powerful  enough  to  be  detectable  at  Earth 
are  generally  large.  Eor  the  15  GLEs  from  cycle  23  for  which  the  area  of  the  parent  sunspot 
group  could  be  determined  (Oliver  2006;  Gopalswamy  et  al.  2012),  the  median  sunspot 
area  of  the  GLE-source  regions  was  750  millionths  of  a  solar  hemisphere  (msh;  range  from 
360-2580  msh,  corrected  for  foreshortening  which  reduces  the  apparent  areas  of  spots  away 
from  disk  center).  As  a  reference,  “naked  eye  sunspots”,  i.e.,  those  that  can  be  seen  without 


Springer 


Solar  Cycle  in  the  Heliosphere  and  Cosmic  Rays 


Fig.  12  A  schematic  view  for 
CME-driven-shock-acceleration 
of  SEPs  indicating  the  locations 
of  quasi-perpendicular  and 
quasi-parallel  CME-driven  shock 
acceleration  of  SEPs  (adapted 
from  Cliver  2009b) 


Quasi-parallel 

Shock 


Ha  Flare 


magnification  through  darkened  glass  or  an  atmosphere  shielded  by  light  clouds  or  smoke, 
have  areas  of  >500  msh  and  make  up  only  about  5  %  of  all  sunspot  groups  (Newton  1958). 

In  the  current  picture  of  particle  acceleration  at  the  Sun  (Reames  1999,  2013;  Tylka  et  al. 
2005;  Tylka  and  Lee  2006;  Cliver  2009a,b),  the  high-energy  protons  that  give  rise  to  GLEs 
are  accelerated  in  solar  flares  and  by  CME-driven- shocks,  either  via  diffusive  acceleration 
at  a  CME  bow  shock  or  via  shock-drift  acceleration  at  the  CME  flanks  as  indicated  in  the 
schematic  in  Eig.  12.  To  drive  the  strong  shocks  associated  with  GLEs,  the  CMEs  need  to  be 
fast  and  massive.  Such  CMEs  dominate  the  energy  budget  of  major  eruptive  flares  (Emslie 
et  al.  2012).  Eor  example,  the  CMEs  associated  with  the  16  GLEs  of  cycle  23  (coronagraph 
data  available  for  15)  had  a  median  speed  of  1810  km/s  (range  938-3242  km/s)  compared 
to  ~420  km/s  for  all  broad  (>60°)  CMEs  in  this  cycle,  while  12  of  the  15  CMEs  were 
full  halos  (Gopalswamy  et  al.  2009,  2012).  The  variation  in  the  occurrence  of  such  CMEs 
(halo  CMEs  with  speeds  >1000  km/s)  for  cycle  23  (1996-2008),  for  which  we  have  the 
homogeneous  SOHO  Large  Angle  Spectroscopic  Coronagraph  (Brueckner  et  al.  1995)  data 
set,  is  shown  in  Eig.  13. 

The  distribution  of  fast  halo  CMEs  in  Eig.  13  looks  like  the  histogram  of  GLEs  over  the 
composite  solar  cycle  for  cycles  18-23  in  Eig.  9,  although  the  broad  maximum  of  the  GLE 
distribution  (spanning  intervals  2-8)  is  shifted  one  year  to  the  left.  The  general  variation  of 
high-energy  SPEs  with  the  solar  cycle  (Eig.  9)  has  also  been  documented  for  significant  SEP 
events  at  lower  energies  (see,  e.g..  Pig.  8).  Comparison  of  Pigs.  13  and  9  (after  normalizing 
for  the  number  of  solar  cycles)  indicates  that  on  average,  only  1  of  about  15  fast  halo  CMEs 
produces  a  GLE.  This  significant  reduction  can  have  multiple  contributing  causes.  Eor  ex¬ 
ample,  increasing  the  speed  threshold  for  Pig.  13  to  1500  km/s,  closer  to  the  ^1800  km/s 
median  speed  for  all  GLE-associated  CMEs  during  cycle  23,  approximately  halves  the  num¬ 
ber  of  fast  halo  CMEs  during  this  period.  Pactors  unfavorable  to  SEP  acceleration  by  fast 
halo  CMEs  include:  low  CME  brightness  (i.e.,  mass/energy)  (Kahler  and  Vourlidas  2005; 
Mewaldt  et  al.  2008),  lack  of  closely-timed  preceding  CMEs  from  the  same  or  nearby  active 
regions  (Gopalswamy  et  al.  2002,  2004;  Kahler  and  Vourlidas  2005),  low  pre-existing  seed 
particle  population  in  the  interplanetary  medium  (Kahler  2001;  Tylka  et  al.  2005;  Cliver 
2006),  and  the  fact  that  GLEs  are  rarely  observed  from  the  eastern  half  of  the  Sun  (only  9 
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Fig.  13  Histogram  of  the 
occurrence  of  fast  (>1000  km/s) 
halo  CMEs  over  the  1 1  equal 
intervals  of  13.3  months  each,  of 
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of  71  cases)  due  to  the  propagation  of  SEPs  along  the  Parker  spiral  (e.g.,  Shea  and  Smart 
1990). 

When  the  parent  solar  activity  (eruptive  flare)  is  near  the  central  meridian  of  the  Sun,  a 
fast  CME  can  propagate  outward  from  the  sun  accelerating  particles  as  the  shock  traverses 
the  interplanetary  medium  between  the  Sun  and  Earth.  When  the  parent  solar  activity  is  lo¬ 
cated  toward  the  western  limb  of  the  Sun  (as  viewed  from  Earth),  the  time/intensity  profiles 
of  these  SPEs  tend  to  have  rapid  rates  of  rise  to  maximum  intensity  followed  by  a  rapid  de¬ 
cay.  Based  upon  those  observations.  Shea  and  Smart  (1994)  and  Smart  et  al.  (2006)  identified 
a  bi-modal  distribution  of  SPEs  classified  as  “interplanetary  shock  dominated  events”  and 
“near-Sun  injection  events”.  Eor  the  former,  a  fast,  broad,  CME-driven  interplanetary  shock 
from  activity  near  the  central  meridian  of  the  Sun  continuously  accelerates  ions  throughout 
its  entire  passage  from  the  Sun  to  Earth.  Eor  the  lower  energy  protons  (e.g.  <  30  MeV)  the 
initial  particle  flux  observed  at  Earth  may  be  relatively  small,  but  the  magnitude  continues  to 
increase  as  the  interplanetary  shock  approaches  Earth,  often  maximizing  as  the  shock  passes 
Earth. 

Eor  the  “near-Sun  injection  events”  there  is  identifiable  solar  activity,  preferentially  on 
the  western  hemisphere  of  the  Sun,  with  an  associated  fast  CME  that  typically  propagates 
through  the  western  heliographic  longitudes  as  viewed  from  Earth.  In  this  type  of  event. 
Earth  is  well  connected  to  the  solar  active  region  by  the  interplanetary  magnetic  field,  and 
the  solar  protons  can  arrive  at  Earth  in  a  prompt  time  frame.  Solar  particle  intensity  observed 
at  Earth  from  activity  on  the  western  solar  hemisphere  typically  has  a  faster  increase  and 
higher  magnitude  than  similar  activity  from  the  central  meridian  or  eastern  sector  of  the  Sun 
(see  Eig.  10  in  Shea  and  Smart  1995). 

4.3  Variation  of  High-Energy  SEP  Events  with  Long-Term  Solar  Activity 

The  era  of  direct  SPE/GLE  recording  is  relatively  short  and  provides  little  information  on  the 
long-term  variability.  On  the  other  hand,  extreme  GEE  events  can  be  identified  as  spikes  in 
production  of  cosmogenic  isotope  (^®Be  and  ^^C)  records  on  the  centennial-millennial  time 
scale  (e.g.,  Usoskin  and  Kovaltsov  2012).  Eor  example,  the  two  largest  SPE/GLE  candidates 
ever  observed  (Miyake  et  al.  2012, 2013;  Usoskin  et  al.  2013;  Oliver  et  al.  2014),  were  found 
in  775  AD  and  993  AD,  as  inferred  from  measured  in  dendrochronologically-dated 
tree  rings.  As  shown  by  Kovaltsov  et  al.  (2014)  such  spikes  are  related  to  the  >200  MeV 
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Fig.  14  Scatter  plot  of  solar 
cycle  totals  of  F200  vs. 

(SSN) cycle  for  cycles  19-23. 
Points  are  also  plotted  for  the  775 
AD  event  and  993  AD  events 
with  the  decadal  SSNs 
encompassing  these  years  taken 
from  Kovaltsov  et  al.  (2014) 


fluence  in  SPEs,  F200,  as  based  on  SEP  spectra  from  Tylka  and  Dietrich  (2009).  Large  F200 
SPEs  indicate  events  with  a  hard  SEP  spectrum.  The  largest  F200  SPEs  since  1950,  i.e., 
those  with  F200  >  3.0  x  10^  proton/cm^,  occurred  on  23-Eeb-1956  (heliolongitude  W80), 
12-NOV-1960  (W04),  15-Nov-1960  (W35),  29-Sep-1989  (^WlOO),  19-Oct-1989  (E09),  and 
14-Jul-2000  (W07).  McCracken  et  al.  (2012)  reported  that  three  of  these  events,  those  that 
originated  at  solar  longitudes  >W35,  exhibited  high-energy  impulsive  peaks  (extending  to 
rigidities  >5  GV)  with  short  (3-5  minutes)  rise  times  and  fall  times.  They  pointed  out  that, 
“Eor  a  quasi-perpendicular  shock  such  as  on  the  flanks  of  a  CME,  it  is  conceivable  that  the 
shock  front  (sweeping  laterally  across  the  face  of  the  Sun)  may  connect  and  disconnect 
with  open  field  lines  [nominally  near  W60]  that  are  connected  to  Earth,  and  so  lead  to 
both  an  abrupt  switch-on  and  switch-off  of  the  SEP  pulse.”  The  three  large  F200  events 
without  a  reported  high-energy  impulsive  peak  all  originated  near  central  meridian.  Eor  such 
events  the  CME-fiank-driven  shock  may  not  reach  the  ~W60  zone  of  good  connection  (or 
may  be  significantly  weaker  when  it  does)  and  the  SEP  event  can  be  dominated  by  the 
outward  expanding  quasi-parallel  bow  shock  of  the  CME.  In  this  case  the  shock  front  can 
remain  connected  for  a  longer  time  to  the  Parker  spiral  to  Earth,  contributing  to  the  extended 
intensity  vs.  time  profiles  observed  for  SEP  events  arising  from  near  central  meridian  (Cane 
et  al.  1988).  The  black-filled  rectangles  in  the  histogram  of  Pig.  9  show  the  distribution  of 
the  six  large  F200  events  over  the  solar  cycle. 

Active  regions  with  sunspot  areas  >1000  msh  can  occur  in  even  relatively  weak  solar 
cycles.  Eor  example,  in  current  low-activity  cycle  24,  the  largest  active  region  yet  observed 
had  an  area  of  1580  msh  (on  4-Peb-2014).  Moreover,  at  any  given  proton  energy,  the  fluence 
from  a  single  large  SPE,  or  a  cluster  of  closely-timed  events  from  the  same  active  region, 
can  dominate  the  total  fluence  for  a  1 1-yr  cycle.  Thus  we  should  not  expect  SEP  fluence  over 
a  cycle  to  be  closely  tied  to  the  strength  of  the  solar  cycle.  This  can  be  seen  for  solar  cycles 
19-23  in  Pig.  14  where  F200  totals  for  all  GLEs  during  a  given  cycle  (Kovaltsov  et  al.  2014) 
are  plotted  against  the  average  daily  sunspot  number  for  that  cycle,  (SSN) cycle-  No  clear 
dependence  of  F200  on  (SSN) cycle  is  apparent  over  the  50-90  range  of  (SSN) cycle  values  for 
these  cycles.  In  the  case  of  the  775  AD  event,  the  F200  fluence  for  a  single  event,  or  cluster 
of  events,  exceeded  the  fluence  from  the  largest  solar  cycle  (No.  19)  during  the  ^1945-1995 
modern  grand  maximum  (Usoskin  et  al.  2007;  Usoskin  2013)  by  about  a  factor  of  thirty. 

Work  by  McCracken  et  al.  (2004)  and  McCracken  (2007b)  implies  that  if  Pig.  14  in¬ 
cluded  smaller  cycles  with  (SSN) cycle  values  in  the  range  from  0-50,  e.g.,  for  the  current 
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cycle  24  (Richardson  2013),  we  might  see  an  anticorrelation  between  F200  and  (SSN) cycle- 
Those  authors  reasoned  that  less  intense  cycles  would  produce  a  weaker  ambient  HMF  and 
suggested  that  CMEs  injected  into  these  weaker  fields  would  produce  shocks  with  higher 
Mach  numbers,  resulting  in  more  frequent  and/or  energetic  SEP  events.  At  the  same  time 
weaker  cycles  should  produce  fewer  energetic  CMEs,  partly  off-setting  the  advantage  of 
lower  ambient  HME  Thus  we  see  a  virtual  absence  of  GLEs  during  the  characteristic  inter¬ 
vals  of  low  HMF  near  sunspot  minima  because  of  the  relative  lack  of  fast  halo  CMEs  at  these 
times  (see  years  1  and  1 1  in  Figs.  9  and  13).  McCracken  and  colleagues  based  their  sugges¬ 
tion  in  part  on  pre-space-age  (e.g.,  cycles  12-14,  1878-1913)  peaks  in  nitrate  concentration 
in  ice  cores  which  they  attributed  to  strong  SEP  events  (McCracken  et  al.  2001,  2004)  and 
also  on  four  intense  GLEs  from  1942-1956  during  a  period  of  low  HMF  strength  inferred 
from  balloon-borne  and  ground-based  ionization  measurements  of  cosmic  rays  (McCracken 
2007a,b).  The  ice  core  evidence,  however,  has  been  disputed  by  Wolff  et  al.  (2012).  In  ad¬ 
dition,  the  low  inferred  HMF  values  circa  1950  inferred  from  ionization  chambers  are  not 
substantiated  by  more  recent  reconstructions  of  HMF  based  on  geomagnetic  observations 
(Svalgaard  and  Cliver  2010;  Lockwood  and  Owens  2011;  Kahler  2008)  as  well  as  on  cos¬ 
mogonic  nuclides  (Steinhilber  et  al.  2010). 

Recently,  the  current  weak  solar  cycle  has  provided  a  more  direct  test  of  the  conjecture 
of  McCracken  and  colleagues — with  mixed  results.  Gopalswamy  et  al.  (2014)  noted  the 
relative  absence  of  GLEs  during  the  rise  phase  of  cycle  24  (see  Table  1).  They  attribute 
this  deficit  to  reduced  total  (magnetic  plus  plasma)  pressure  in  the  heliosphere  that  leads, 
in  turn,  to:  increased  CME  lateral  expansion,  dilution  of  CME  magnetic  fields,  and  less- 
efficient  shock  acceleration  of  SEPs.  However,  Gopalswamy  et  al.  (2014)  also  note  that 
while  the  number  of  GLEs  thus  far  in  cycle  24  is  significantly  reduced,  the  overall  number 
(28)  of  significant  >10  MeV  proton  events  is  comparable  to  that  (34)  for  the  corresponding 
phase  of  the  stronger  cycle  23.  Following  McCracken  and  colleagues,  they  suggested  that 
the  similarity  in  the  counts  of  the  lower  energy  SEP  events  for  the  two  cycles  may  result 
from  the  relative  ease  of  shock  formation  due  to  the  lower  ambient  solar  wind  Alfven  speed 
in  cycle  24,  off-setting  the  size  advantage  of  cycle  23. 

The  current  low  cycle  is  providing  valuable  data  and  insight  on  SEP  acceleration  during 
conditions  of  low  solar  activity.  While  significantly  more  active  than  the  quietest  ~50  years 
of  the  Maunder  Minimum  (e.g.,  Eddy  1976;  Usoskin  2013),  cycle  24  is  giving  us  our  best 
glimpse  to  date  of  SEP  activity  during  periods  of  low  solar  activity. 


5  Summary 

In  this  concise  review  we  discuss  the  solar  cycle  and  its  manifestations  in  the  heliosphere 
and  cosmic  rays. 

Most  heliospheric  parameters  are  measured  in  situ  near  Earth,  and  thus  do  not  straight¬ 
forwardly  represent  the  entire  heliosphere.  While  some  parameters  (e.g.,  HMF)  are  roughly 
representative,  others  (solar  wind  velocity  or  plasma  density)  are  not  representative  for  the 
whole  heliosphere.  Global  heliospheric  indices,  such  as  the  open  magnetic  flux  are  useful  in 
this  respect.  Galactic  cosmic  rays  can  serve  as  probes  of  the  3D  heliosphere,  reflecting  its 
variability  on  the  time  scale  of  11 -year  cycle  and  beyond.  The  quasi  22-year  variability  in 
GCR  observed  as  alternation  of  top-  and  flat-peaked  cycles  is  caused  by  the  drift  effect  of 
GCR  transport  in  the  heliosphere.  Special  attention  is  paid  upon  the  recent  minimum  of  solar 
activity  with  very  quiet  heliospheric  conditions  (low  values  of  HMF  and  the  interplanetary 
turbulence)  which  led,  together  with  the  flat  HCS,  to  the  record  high  level  of  GCRs  near 
Earth  in  2009. 
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We  also  discuss  the  statistics  of  the  occurrence  of  SPEs  and  GLEs  that  serve  as  a  probe 
for  the  inner  heliosphere  and  solar  coronal  activity.  The  occurrence  of  major  events  depicts 
an  overall  tendency  to  follow  the  solar  cycle  but  individual  events  may  appear  at  different 
phases  of  the  solar  cycle,  since  it  is  defined  not  only  by  the  solar  energy  releases  but  also 
by  the  location  of  these  releases  relative  to  the  nominal  ~W50°  footpoint  of  the  magnetic 
fieldline  that  connects  to  Earth.  We  also  discuss  that  the  occurrence  of  major  SPEs  is  not 
directly  related  to  the  overall  level  of  solar  activity,  and  strong  events  may  occur  even  during 
moderate  solar  cycles. 

In  summary,  the  solar  cycle  in  the  heliosphere  and  cosmic  rays  depicts  a  complex  pattern 
which  includes  different  processes  and  cannot  be  described  by  a  simple  correlation  with 
sunspot  number. 
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